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Abstract 

Cosmic strings arising from GUTs can catalyse baryon decay processes 
with strong interaction cross sections. We examine the mechanism by 
which the cross section is enhanced and find that it depends strongly on 
the details of the distribution of gauge fields within the string core. We 
propose a calculational scheme for estimating wavefunction amplification 
factors and also a physical understanding of the nature of the enhancement 
process. 
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1 Introduction 



Grand unified theories (GUTs) predict a rich variety of topologically stable "de- 
fects" - domain walls, monopoles, cosmic strings - whose localized concentrations 
of unbroken gauge fields and Higgs condensate would be expected to catalyse 
baryon decay. In particular, there would be a non-zero amplitude for quarks 
which penetrate the core of such a defect to decay into leptons, as quarks and 
leptons appear in the same multiplet in GUTs. The typical size of such de- 
fects (monopole radius, or cosmic string radius) will be of the order of the GUT 
length scale (10~^''cm). Naively, one might expect the cross section for monopole 
catalysed baryon decay processes to be of the order of the corresponding area 
(10~^''cm^) and that for decay processes mediated by a string, the cross section 
per unit length of string would be of the order of the GUT length. In the case of 
monopoles however, it has long been understood, since the work of Callan and 
Rubakov [|l|, |^, that cross sections for such processes will be greatly enhanced 
and actually be of the order of the QCD area (10"'^° cm^). The essential point 
is that there is a mechanism (due to the long range external magnetic field of 
the monopole) for the amplification of the quark wavefunction at the monopole 
core and thereby an increased probability of penetration. Since quark masses 
and energies will be of the order of the QCD scale which is typically 15 orders of 
magnitude smaller than the inverse radius of the monopole, one is thus interested 
in the low energy scattering of quarks on monopole targets. (We shall, however, 
treat energies sufficiently high that quarks may be treated as free.) 

In the case of cosmic strings, mechanisms for a similar enhancement of baryon 
catalysis have often been discussed 0,^0. Afford and Wilczek showed that 
a GUT string can carry fractional flux in units of 27i/q, where q is the charge of 
a quark or lepton in the theory. Such a particle thus has a significant low energy 
elastic cross section for scattering from the string due to the presence of the 
topologically non-trivial, but pure gauge external field configuration, as shown 
long ago by Aharonov and Bohm 0. It was also noted in that the fermion 
wavefunction is amplified at the string. Physically, this means that the presence 
of the external gauge fields has increased the probability of fermions penetrating 
the region of unbroken symmetry, thereby enhancing the catalysis rate. 

In this paper, we follow the computational scheme for catalysis cross sections 
originally proposed in ^ which breaks the calculation into two steps. In the 
first, the decay cross section is computed using free fermions (i.e. ignoring the 
external gauge fields) to give a GUT scale cross section, known as the geometric 
cross section. This step clearly depends strongly on the decay mechanism - 
whether baryon decay is mediated by interactions with internal X and Y gauge 
fields or by interaction with a scalar condensate in the string core. This is partly 
determined by the model of the string used. In the second step, one computes 
the degree to which the geometric cross section is enhanced by considering the 
scattering of fermions off the string. The prescription used is derived from first 
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order perturbation theory and gives the decay cross section as 



da ^4 da 

dn~ dn 



(1) 



geom 



where {da/dfl)geom is the geometric cross section and A is the wavefunction am- 
plification factor 

^ |^(a)| . . 

|^free(a)| ^ ^ 

i.e. the ratio of the magnitude of the spinor in the magnetic field to the free 
spinor evaluated at a the core radius of the string (at infinity, both spinors are 
normalised by scattering boundary conditions). 

It was shown in 0] that the total cross section can be enhanced up to the 
QCD scale. The second step also depends on the internal model of the string 
(gauge fields or scalar condensate) and also on the net flux carried by the string; 
however, it was concluded in that the distribution of gauge fields within the 
core does not affect the amplification factor. This claim was made on the basis of 
a consideration of two simple models of flux distribution: the case of a flux ring 
at the core radius and the case of uniformly distributed flux within the core. 

Here, we re-examine this claim by considering more general models of the flux 
distribution. We shall concentrate exclusively on the case where the dominant 
low energy scattering of fermions off the string is due to interactions with gauge 
fields and where the decay process itself is mediated by gauge fields, although our 
general methodology could easily be extended to cover interactions with scalar 
fields. Under certain assumptions, we demonstrate that the results of are 
indeed independent of the details of the flux distribution. However, when the flux 
is allowed e.g. to change sign in the string interior, we find that the amplification 
factors can be strongly dependent on the details of the flux distribution. It might 
be objected that such a field configuration is unphysical, and certainly a single 
gauge field whose flux changed sign in the core would probably be unstable. 
However, the situation we envisage is where two or more gauge fields may be 
represented by a single effective U{1) gauge field. For example, these fields might 
be the X and Y gauge fields and the electromagnetic field. If the separate gauge 
fields have different ranges (as they do in this example) and if they are of opposite, 
but constant sign (and there is no a priori reason to prevent this) then the 
effective gauge field could certainly change sign in the core without prejudicing 
the stability of the string. Thus our results may be important for computations 
with realistic string models. 

In order to treat these problems, we develop and extend a calculational scheme 
(the scattering length formalism) which was originally developed |^ (see also 
|T^ ) to study the general problem of the large scale behaviour of small objects. 
A remarkable property of Dirac operators coupled to external U{1) gauge fields 
allows us to calculate the relevant parameters (the scattering lengths introduced 
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below) analytically for arbitrary flux distributions and by developing the analogue 
to the low energy expansion of potential scattering theory [|12| we are thus able to 



provide a simple means of estimating the ampliflcation factors for baryon decay. 
Our formalism leads to a physical interpretation of the enhancement process. 
We also comment on the validity of the thin wire approximation to the flux 
distribution used in 0, Q and relate our results to our other work on the large 
scale effects of small objects P, ||. 

Our principal assumptions are as follows: 

1. Quarks are treated (as in [|, |^) as free Dirac particles with energies above 
the conflnement scale. (On the GUT scale, this still corresponds to low 
energies.) 

2. Decay processes are mediated by interactions with gauge flelds in the core. 

3. The cosmic string is assumed to be an infinitely long straight cylindrical 
string along the 2;-axis of radius a. 

4. All fields are cylindrically symmetric about the z-axis. 

5. The quark wavefunction is z-translationally invariant. (Due to the low 
energies of incoming quarks, we expect that including the z-dependence of 
the quark wavefunction will not significantly alter the physics.) 

Our conventions are as follows: the metric has signature H , the in- 
coming quark has charge — e and the electromagnetic vector potential is defined 
by = (0, —A) with V A A = B, where B is the magnetic fieldQ. (r, 6, z) are 
cylindrical polar coordinates about the 2;-axis. 



2 The Dirac Equation with a Flux Tube 

We consider the minimally coupled Dirac equation (z7^(9^ — ieA^) — m)il}i = 
in the 7-matrix representation 

For a fiux tube of magnetic fiux B = B{r)z with B{r) vanishing for r > a, a 
simple Stokes' theorem argument gives 

A = ^e 

er 

^ We refer to in terms appropriate to an electromagnetic field; however it should of course 
be thought of as the effective gauge field. 
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where a{r) is defined by 



a(r) 



B(r')r'dr'. 



We define $ = a (a) = (total magnetic flux)/ (27r/e). In addition, we define 
= $ — [$], where [$] is the greatest integer strictly less than $. We will only 
be interested in the case of non-integer flux. 

Diagonalising Jz, the angular momentum operator about the 2;-axis, given by 



J. 



1 / 0^3 



(3) 



with eigenvalues n + j {n El) using the ansatz 



ijji = e exp < i 



^1 1 / ag 
"+2-2 aa 



(4) 



we find that the Dirac equation separates into two 2-spinor equations for Fj, G\ 
and FJ^, with radial equations 



-i[uo — m 



(5) 



where m = m for Fj , Gj^ and m 



-m for F^,^^,. 



In the rest frame of the particle, t and | correspond to spin aligned and anti- 
aligned with z. We now drop the arrows and proceed to treat only one 2-spinor 
in each angular momentum sector: 



Fn{r) 
GJr) 



(6) 



It will turn out that the amplification factor has the same order of magnitude 
whichever 2-spinor is chosen. 

The equations decouple to give 



1 d d 2 
— rf^-r + Vn +"^1 
r dr dr 



I d d 2 21^ 

-—r— + V-+m'-u^\Gn 
r dr dr 








where the effective potentials are given by 



n + a;(r)\ a'{r) 



+ 



(7) 
(8) 

(9) 



In the subsequent discussion we shall see that an understanding of the mech- 
anism for enhancement depends crucially on a consideration of the details of the 
effective potentials (i.e. on the details of the flux distribution) inside the 
cosmic string, and that one loses essential insight if one passes to the thin wire 
approximation, or if one takes non-smooth models for the flux distribution. 

Equations (|^) and (|D are not independent, but are coupled via the first order 
equations (|]). It is easy to see that the coupling ensures that if one component 
is regular at the origin, then so is the other. Solving the system at momentum 
k = \JuP' — fh? with regular boundary conditions at r = 0, we find that outside 
the string core, where a{r) = <I>, -F„ ^(r) and Gn^ki.'f) may be written as 



where On and must be determined by matching the external solution to the 
solution inside the string core. The coupling between equations (||) implies that 
ianOn = —iasiipn in all sectors except n = — 1 — [$], where we require tan^„ = 

-cot ipn- 

If we were to employ the "thin wire approximation", in which the flux is 
concentrated in a flux tube of infinitesimal radius, the external solutions ([TT|) , (p!^) 
would hold down to r = 0. We could then apply the criterion of local square 
integrability (with measure rdr) at r = to fix 6'„ = = (for all k) in all 
sectors other than n = —! — [$], which we refer to as the critical sector. In 
the critical sector any choice of 9n and ipn consistent with the coupling tan^„ = 
— cot ipn leads to a locally square integrable wavefunction. We therefore see that 
in the thin wire approximation, all wavefunctions are regular and vanishing at the 
origin except in the critical sector where there is a 1-parameter family of possible 
boundary conditions, for each of which at least one of the 2-spinor components 
must be irregular at r = 0. 

Mathematically, this corresponds to the fact that the Hamiltonians Hl^ de- 
rived from (^) are essentially self-adjointQ on C°°-spinors compactly sup- 
ported away from r = in L'^{{0, oo)'^, rdr) in all sectors other than the critical 
sector where Hl^ have deficiency indices (1, 1) 0, This entails that there is a 
1-parameter family of self-adjoint extensions labelled by elements of U{1), each 
of which corresponds to a different choice of boundary condition in the critical 
sector. The self-adjointness of the Hamiltonian is necessary to ensure a unitary 
time evolution, so it is only by choosing a particular self-adjoint extension that 
we can specify a well-defined global dynamics for the system (see also [§, H). 

Following the procedure for computing amplification factors introduced in 
(see equation (^), we now impose scattering boundary conditions on the 
wavefunction in each sector and then compare the magnitude of the wavefunction 

^An operator A on some domain is essentially self-adjoint if its operator closure A is self- 
adjoint i.e. A and its adjoint A* have the same domain on which they act in the same way. 



Fn,k{r) oc cos6'„J|„+$|(fcr) - sin 6'^ J_|„+$|(fcr) 
G„,fe(r) oc cosLpnJ\n+'S>+i\ikr) - sinv9„J_|„+$+i|(A;r) 



(11) 
(12) 
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ai r — a against that of the free wavefunction, which is of order 1. The scattering 
boundary conditions are derived in Appendix A and give the normahsed spinor 
as: 



1 _ (-i)2|n+*ltan^„] (-i)l"+*lj|„+$|(fcr) 

1 - (-z)2|"+*+^ltan(^J"' A(-2)I"+*+i|J|„+<,+i|(A;t) 



+ 



1- f-i)-2|"+*lcot^. 



i+<i>| J 



|n+$| 



{kr) 



1 _ (-i)-2|n+*+ilcot(^J" A(-i)-l-+*+i|j_|„+$+i|(A;r) 



(13) 



where A = —k/{u} + m) which is of order 1 at the energies of interest (where 
k m), for either choice m — ±m. 



3 Scattering Length Formalism 

In order to determine the wavefunction amphfication, it now suffices to specify 
6n and This is accomphshed by performing the analogue of the low energy 
expansion in potential scattering theory. If we denote the logarithmic derivative 
K,k/Fn,k\r=a by -Dn, wc may expand £>+ in powers of {kaf: D+ = + 
{kayD^^'^ + 0{ka)^. Note that ka, the product of fermion momentum and string 
radius is of order 10~^^. The matching between internal and external solutions 



cot 9n 



\n+^\{kr) </i|„_|_$|(A;r) 
D„J|„+$l(A;r) - J|;+$|(fcr) 



(14) 



(where the prime denotes differentiation with respect to r) may be expanded to 
give 



cot 6r, 



r(i + |n + $|) 



l-^-^^^{kr:r + 0{kA (15) 



^ 2 J r(l - |n + $|) [ 2 

where we call the scattering length of the effective potential given by 

"aDlo)+ - In + $ 



Rt 



l/{2|n+$|) 



(16) 



.aL'l°^++ |n + $, 

and r^, which generalises the effective range of potential scattering, is given by 



' n I 



2'K\n + $1 



-2|n+*| 



{ (^) 



\n + <l>| 



2 + 



2|n+*| 



1 + |n + $| 



/ a 



— |n+$r 



(17) 
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Our parameter 6n is related to the phase shift Sn relative to the Aharonov-Bohm 
scattering (defined by oc cos5„J|„+$|(/i;r) — sin (5„A^|„4.$|(/i;r)) by 



cot 6r, 



cos(|n + $|7r) — cot 9n 
sin(|n + $|7r) 



Note that there is another definition of scattering length (which is the one 
used in p, ^): equation (|^) may be solved exactly at a;^ = for r > a (where 
Vn{r) = (n + $)^/r^) and takes the simple form 



Fn(r) oc 



(19) 



which we can use to define R^. Note that (to ensure reality of F„(r) up to a 
phase) (i?^)^'""'"*' must be real (although possibly negative) and so the allowed 
values of R'^ lie on a contour in the complex plane. We make this definition of 
R^ to ensure that our 'scattering length' really has dimensions of length; it will 
turn out that this is a natural parametrisation to use. i?^ may be expressed in 
terms of Fn{r) by means of a fitting formula 



R1 



all- 



2\n + <t>\ rl"+*lF„ 



(rl"+*lFO' 



l/(2|n+$l) 



(20) 



Substituting = F'^{a)/Fn (a) in (|20| ) reduces it to (|16D and we see that our 

two definitions of scattering length agree. 

Similar considerations for cot (when the logarithmic derivative = 
G'^/Gn\r=a is expanded as = D^^)" + {kafD^:^^- + O(fca)^) give 



cot ipn 



'kR- 



-2|n.+<I>+l| 



r(l + In + 'I' + ll) 
r(l-|ra + <l> + l|) 



1 - 



In + $ + llvr 



{kr-f + 0{kY 



(21) 

where i?„ is determined from the zero energy external solution Gn{r) by the 
fitting formula 



R- = a^l 
and r~ is given by 



2|n + <l> + l| rl"+*+ilG'„ 



27r|n + $ + l| 



{a/R; 



-2|n+$+l| 



l/(2lr!,+'I>+l|) 



(22) 



(«/^n) 



2|n+$+l| 



+2aD, 



(1)- 



a 

rE 



1 - |n + $ + 1| 

\ In+'I'+ll 



1 + In + $ + 11 



a 
R^ 



(23) 
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In this case, the zero energy external solution is 



^ \ |n+*+l| / „ \ -In+'I'+ll 




Gn{r)^\ — ] -( — ) (24) 



The remarkable feature of equations (^ and (P) to which we referred in the 
introduction is that they may be solved analytically at zero kinetic energy (cj^ = 
iv?). This is because they factorise at zero kinetic energy (see Appendix B) 
and may be viewed as a consequence of an abstract supersymmetry possessed by 
Dirac operators coupled to external magnetic fields (see |]TT|). In Appendix B, we 



derive these solutions and place bounds on their corresponding scattering lengths 
sector by sector. We summarise our results in Table 1. In each sector, one of 
the two scattering lengths is either zero or infinite and so the corresponding 
low energy expansion ( ([T5|) or (pT]) ) breaks down. One could derive the form of 
the expansion in these special cases; however, the other low energy expansion 
remains well-defined and so, by using the formulae connecting On and we can 
always determine both as functions of k at low energies. 

In sectors where the bounds derived in Appendix B allow a range of possible 
scattering lengths, the precise details of the flux distribution fix a particular 
choice via the fitting formulae. The bounds derived in Appendix B can be shown 
to be 'best possible' and so, in sectors where our bounds permit, can be made 
arbitrarily large for arbitrarily small values of a. This persists in the limit as 
a — > contradicting the results of IT^ . 

Large scattering lengths occur only when the corresponding effective potential 
exhibits a potential well, as may be seen by the following argument. If the effec- 
tive potential in, say, equation (^ is everywhere non-negative, then a convexity 
argument applied to the differential equation (see 0) shows that the zero energy 
solution Fn{r) ^ with regular boundary conditions at the origin must satisfy 
F^{a)/Fn{a) > 0. Inserting this in the appropriate fitting formula (^) we find 
that the corresponding scattering length satisfies < a. Thus to generate 
large scattering lengths it is necessary for the effective potential to exhibit a well. 
Note that if $ > 0, it is necessary that the magnetic field B{r) change sign within 
the core of the string for i?+ to be large in the critical sector; for the existence 
of a well in implies that a'{r)/r = B{r) must be negative in some interval 
in (0,a). Similarly, for $ < 0, R~ is large only if B{r) changes sign. We have 
already observed that such configurations are not necessarily unstable when an 
effective gauge potential is considered. 

Of particular interest in Table 1 are the cases —[^]<n<—l for $ > 1 and 
< —2 — [$] for $ < — 1. The infinite scattering lengths in these sectors are due to 
the presence of bound states of zero kinetic energy located at = — e($)m (where 
e(x) = ±1 cLS X IS greater than or less than 0). In accordance with a theorem 
of Aharonov and Casher [|T0|, there are precisely [$] such states for $ > and 
— [$] — 1 such if $ < 0. Although one of the scattering lengths in the critical 
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sector n = — 1 — [$] is infinite, this is not a bound state, as the wavefunction at 
zero kinetic energy fails to be square integrable at infinity. (RecaU: [$] is the 
greatest integer strictly less than $.) 



4 Calculation of Amplification Factors 

The information in Table 1 and expansions (|T5|) and (|2T|) enables us to gain 
some insight into the form of the wavefunction amplification at low energies (in 
comparison with the GUT scale), as we now have some control over the leading 
order behaviour of cot and cot<^n. We can use this to construct simple order 
of magnitude arguments which are sufficient to demonstrate the range of possible 
behaviour. In particular, we will assume that aDj^^"^ have magnitude of order 
1 or smaller - we will not consider the additional behaviour occurring if -D^^-*^ 
are finely tuned so as to produce cancellations in (|1^ and (pSj). We also ignore 
the effect of higher terms in the low energy expansions (p!5D and (|2l|) . These 
assumptions amount to the approximation 



O ((A;a)2(a//2+)2|«+*l] i?+ < a 
{kr^r^{ 0{{kaf) Rt-a . (25) 

o((A:a)'(^;[/«)'i"+*i) ^;[»« 

This allows us to conclude that cot 6 is well-approximated by the first term in (p!5|) 
provided that 

a(A;a)^/l"+*l < < a(A;a)-i/l"+*l (26) 

and so cot 6'„ = 0((A;-R^)~^l"+*l). In the language of we say that is 
'believable' at scale k~^. For very small scattering lengths R^ <^ a(A;a)^/l"+*l, 
we have cot0„ ^ (fca)^^*-"'^^'"^*'^ and for very large scattering lengths R^ ^ 
a(A;a)-^/l''+*l, we have cot 6* ~ (A;a)2(i-l"+*l). Note in particular that QCD scale 
scattering lengths i?^ = 0{k~^) are classified as 'very large' unless |n + $| < 1 
when they fall within the range (pB]). Similar conclusions may be derived for 
cot ip. 

We treat the case ^ > and distinguish two cases: the critical sector n = 
— 1 — [$] and other sectors with n 7^ — 1 — [$]. In the critical sector, we find 
(inserting the coupling relation tan On = — cot ifn in (|T3|)) 



+ 



. ( . 1-1/ HV-'Ji-.ikr) \ 



and so we see that the wavefunction is always amplified at r = a, regardless of the 
value of 6n- However, the degree to which it is amplified is determined by 6'„ and 
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can vary between (ka)"" and {ka)""^, depending on the details of the internal 
flux distribution. As |n+ $| < 1 in this sector, we find that we can approximate 
the expression (|1^) for cot6'n by its first term over a large range of values for 

including QCD length scales. We find that for 'very small' scattering lengths 
or values of within a few orders of magnitude of the tube radius a (i.e. for 
GUT scale scattering lengths) the amplification is of order {ka)~^ and will occur 
in the lower component of the spinor. For scattering lengths on the QCD scale, 
for which fci?^ ~ 1, the amplification is given by the larger of {ka)~^ (lower 
component) and {kaY~^ (upper component), whilst for i?+ in excess of the QCD 
scale, the amplification is of order {kaY~^, occurring in the upper component. 
Thus the amplification factor is strongly dependent on the details of the flux 
distribution. 

In sectors other than n = — 1 — [$], we find 



1 - (-z)2|"+*ltanft 
1 - (-z)2|"+*l cot 9. 



-1 / (-^)l"+*lj|n+*|(A:r) \ 
1^ A(-Ol"+*+V|„,+$+i|(fcr) I 

(-0-l"^*l^-|n+*|(^r) \ 
1 A(-z)-l"+*+i|j_,„^^+.(A:r) i' ^"^^ 



We treat the case — [$] < n < —1 first which arises only for $ > 1 and 
corresponds to the Aharonov-Casher states. Here, we have < i?^ < a. For 
cases in which i?+ is of the order of a, we see that cot^„ = 0((A;a)~^'""''*') as 
{kr^Y <^ 1 and so 



^ , 0{{kat^''\) \ _ ( 0((HI"^*I; I .291 
I 0((A;a)""'"*'"'"+*"'") ; V 0((fca)""'*~'^ ' ^ ' 



and we therefore find amplification only for n = — [$], where the amplification 
factor is (ka)""^ in the lower component. However, if is 'very small' so that 
(kr^Y is no longer negligible, cot^^n ^ (A;a)^^*^-'^+l"+*l) and it is easy to see that 
there is no amplification even in the sector n = — [$]. 

We now treat the case n < — 2 — [$] . In these sectors it is easy to see, using the 
same methods as above, that 'very small' scattering lengths or scattering lengths 
within a few orders of magnitude of the core radius a give no amplification. How- 
ever, if R^ is of the order of the QCD scale or larger, cotft„ = 0((A;a)^(^~l'^"^*l)), 
giving 

0((A;a)l"+*l-2l 
0{{ka) 

and so there is amplification only in the sector n = —2 — [$], with amplification 
factor {ka)~^ in the upper component. The possibility of amplification from this 
sector does not appear to have been noted before. To conclude the analysis for 
$ > 0, it remains to consider n > 0. In these sectors, i?^ is forced to be zero 
and we must therefore consider R~ to give cot ipn and then use cot 6n = — cot ipn- 
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As before, we find no amplification if is 'very small' or within a few orders of 
magnitude of a. In the case where R~ is QCD scale or larger, we find cot (fn = 



Thus there is no amplification for n > unless [$] = in which case, there is 
amplification of order {kay~^ in the lower component in sector n = only. Note 
that in this case, there are no Aharonov-Casher states, and that = is adjacent 
to the critical sector. 

We can derive the analogous results for $ < by sending $ — ^ — $, z/ — >■ 1 — j/, 
n — ^ — 1 — n and R'^ R~ . We find that the critical sector n = — 1 — [$] always 
exhibits amplification: for values of R~ within a few orders of magnitude of a or 
smaller, amplification occurs in the upper component with factor {kaY^^, whilst 
for R~ in excess of the QCD scale, the lower component is amplified by factor 
(A;a)"^. In the Aharonov-Casher states < n < —2 — [$], we find that only 
n = —2 — [$] contributes with amplification (ka)'" in the upper component 
unless -R~ <^ a(/ca)^/'"'"'"*"''^' when there is no amplification. The sector n = — [$] 
allows amplification {ka)"^^ of the lower component only for large R~ and there 
can be amplification in n = —1 only for [$] = —1, when the upper component is 
amplified by {ka)~^ only for large i?^. No other sectors contribute. 

In conclusion, and subject always to the provisos stated before equation (^) 
there are thus at most three sectors which can contribute to wavefunction ampli- 
fication: the critical sector n = — 1 — [$] (in which there is always amplification) 
and the two adjacent sectors. We have also seen that, of the sectors adjacent 
to the critical sector, only the Aharonov-Casher sector (when present) ampli- 
fies generically and requires an anomalously small scattering length to be sup- 
pressed, while non- Aharonov-Casher sectors require anomalously large scattering 
lengths in order to contribute. We summarise our results in Figure 1, where we 
graph p($), which determines the overall amplification factor of equation (^ as 
A = (ka)~^^'^\ Figure 1(a) shows the maximum (solid line) and minimum (dot- 
ted line) possible amplification for each flux. In the most general case, when we 
make no assumptions about the form of the flux distribution, we can say no more 
than this without explicitly computing the relevant scattering lengths. 

However, if we know that the flux distribution is single-signed within the 
core, then if $ > we have V^^ > for all n and so our earlier arguments show 
that < R^ < a; conversely, if $ < 0, we know that < R~ < a. In either 
case, (and provided |$| > 1) the relevant scattering lengths in the contributing 
sectors are bounded between and a. (In the case |$| < 1 (in which there are no 
Aharonov-Casher states), there is a possible contribution in the n = [n = —1) 
sector for $ > ($ < 0) which is due to Rq (R-i) and therefore unaffected 
by these bounds.) Thus for a sufficiently 'nice' subclass of single-signed flux 



0((fca)2(i-l"+*+i|)) and so 




(31) 
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distributions, the scattering lengths of interest will all be within a few orders of 
magnitude of a (i.e. of the GUT scale) and so follows the graph shown 
in Fig. 1(b) as only the critical sector and (if |$| > 1) the adjacent Aharonov- 
Casher sector contribute. The solid lines indicate the range of $ for which the 
critical sector provides the dominant contribution to the amplification, whilst the 
dotted portions indicate the ranges where the adjacent Aharonov-Casher state 
gives the dominant amplification. This situation holds for many simple models 
of flux distribution (in particular for those examined in |P) and, as we have 
indicated, for all sufficiently 'nice' single-signed flux distributions. The results of 
correspond to (and agree with) our results in this case and so our discussion 
has demonstrated the extent to which those results can be considered generic. 

We note that for |$| > 1/2 Fig. 1(b) follows the maximum amplification plot 
and so the effect of abnormally large or small scattering lengths could change 
the amplification factor only by suppressing it, as would occur, for example, if 
the scattering length in the critical sector was large (QCD scale or larger), or 
the scattering length in the adjacent Aharonov-Casher sector was much smaller 
than a. When |$| < i, however, it is possible to increase the amplification factor 
considerably by tuning the scattering length in the critical sector to be large (a 
similar effect occurs in this case if the scattering length in n = ($ > 0) or 
n = — 1 ($ < 0) is tuned to be large). The special status of |$| < ^ in Fig. 1(b) 
is due to the absence of an Aharonov-Casher state in this case, which provides 
the dominant amplification when |$| > 1 and the fractional part of |$| is less 
than |. 

Finally, Fig. 1(c) shows the results which would be obtained using the thin- 
wire approximation, on the assumption that the scattering length in the critical 
sector (which is the only free parameter) is of the GUT scale. The inadequacy 
of this approximation is seen by its disimilarity to Fig. 1(b) and the importance 
of the adjacent Aharonov-Casher state becomes clear. It is important to note 
that the thin-wire does not support Aharonov-Casher states as they fail to be 
normalisable. 

5 Conclusion 

We first consider the relation of our current results to our other work p, ^ 
on the large scale effects of small objects. In P] we point out that in many 
physical situations, a small object may be replaced by a point-like or line-like 
idealisation and that if the dynamics of the idealised system admits more than 
one consistent choice of boundary condition (in our case, the Hamiltonian fails 
to be essentially self-adjoint on a suitable domain), this is often a signal that the 
large scale behaviour may be sensitive to the details of the internal structure of 
our original small object. Furthermore, in such cases, the large scale dynamics 
of the true system is well-approximated by the idealisation with an appropriate 
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choice of boundary condition (here a self-adjoint extension) and that therefore 
the parameter(s) labelhng the choice of boundary condition (in our case, the 
scattering lengths) parametrise the possible large scale behaviour. This is the 
content of the "principle of sensitivity" enunciated in pj. 

In the case at hand for the specific purpose of computing wavefunction am- 
plification factors the thin wire approximation fails to be a good idealisation of 
the true system because there are contributions from non-critical sectors, gener- 
ically from the adjacent Aharonov-Casher sector. We note however, that the 
Aharonov-Casher states appear for a quite special and deep reason: an index 
theorem related to the abstract supersymmetry of the Dirac operator, and so 
confirms our general philosophy in that when the principle of sensitivity fails 
to apply, it fails for 'interesting reasons'. 

However, we see that the scattering length formalism developed in Q is 
still applicable and that the amplification factor in the critical sector is strongly 
sensitive to the details of the internal flux distribution. Moreover, if one considers 
the elastic scattering cross section rather than wavefunction amplification factors, 
it is found that the main deviation from the pure Aharonov-Bohm cross section 
at low energies (large scales) occurs in the critical sector and is parametrised by 
the scattering length there. Also, one can show that if a sequence of Dirac 
Hamiltonians describing flux tubes of steadily decreasing radius has a limit (in a 
suitable sense of convergent dynamics, technically strong resolvent convergence) 
which is self-adjoint (i.e. a well-defined limiting dynamics) then the limit must 
be a self-adjoint extension of the idealised thin-wire approximation. Thus the 
principle of sensitivity seems to apply as far as scattering cross sections are con- 
cerned. 

We also note that if one modifies the Hilbert space or the domain on which the 
Hamiltonian is defined in an appropriate way, it is possible to arrange that the 
sectors n = — 2 — [$], n = —! — [$] and « = are precisely the sectors in which 
the Hamiltonian fails to be essentially self-adjoint and that therefore it might be 
that our results can be reconciled with a discussion of self-adjoint extensions after 
all. This may be done in a variety of ways; for example by taking the Hilbert 
space to be the Sobolev space given by the completion of the space of smooth 
spinors compactly supported away from the flux line in the norm defined by 
{(p \ ip) = {(J) \ H'^ijj)L2, where H is the thin wire Hamiltonian. Alternatively, 
one can keep the original Hilbert space, whilst restricting the domain of the 
Hamiltonian to be the range of the massless thin wire Hamiltonian acting on 
smooth spinors compactly supported away from the flux line. We hope to return 
to the significance of these modified versions of the thin wire approximation 
elsewhere. 

We conclude with various remarks. Firstly, the above arguments have estab- 
lished that the amplification factors for baryon decay enhancement calculations 
can depend substantially on the internal distribution of the magnetic flux. In 
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particular, we note that the case |$| < ^ in which we have seen that amphfica- 
tion can be increased includes two of the most physically interesting cases in the 
GUT model of 0, where scattering of the d quark is modelled by $ = — |, and 
the electron as $ = ;|. For these values of $, the difference between {ka)~^ and 
{kaY~^ amounts to 8 orders of magnitude. It is therefore of some importance 
that the amplification factor be correctly computed, taking into account the de- 
tails of the model. The scattering length formalism presented here provides a 
convenient calculational scheme. 

Secondly, we turn to the physical interpretation of the process of wavefunc- 
tion amplification. From above, it is clear that for a particular component to be 
amplified it is necessary that its corresponding (zero energy) scattering length be 
large (QCD scale or above). (That it is not sufficient may be seen by considering 
the Aharonov-Casher states). We also saw above that large scattering lengths 
occur only when the effective potential exhibits a well of negative potential. This 
makes it reasonable to suggest that the physical cause of wavefunction amplifi- 
cation (and therefore of baryon decay enhancement) is a resonance phenomenon 
caused by the spin-fiux interaction: incoming quarks may tunnel into the well 
and be delayed, perhaps being refiected by the walls of the well before tunnelling 
out. Quarks are therefore present in the core of the string much longer than 
would naively be expected and therefore decay processes occur with increased 
probability. 

This interpretation of the enhancement process as a resonance phenomenon 
depends on an examination of the details of the effective potentials (^ and (^) 
- in particular the presence or absence of wells. Thus our interpretation did not 
emerge clearly from previous work on this subject, where the restricted range of 
particular models treated did not display all of the possible qualitative features 
discussed above. We have seen in particular that the thin wire approximation 
(and therefore an approach based solely on self-adjoint extensions on the usual 
domain) is inadequate for this problem, as we have found possible contributions 
to enhancement not only from the critical sector n = — 1 — [$] (which provides 
the only contribution in the thin-wire approximation) but from the two sectors 
adjacent to this sector. In particular the adjacent Aharonov-Casher sector pro- 
vides the dominant enhancement for certain ranges of $ (provided the associated 
scattering length is of the order of a, which is the case e.g. for the simple fiux 
distributions models of |0). This relation with the Aharonov-Casher state has 
not been noted before. 

Finally, we note that although we find one more contributing sector than ||^, 
the amplification arising from the new sector is at most of the order of that from 
the other two sectors. This is in accord with the unitarity bounds established in 

i- 
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A Scattering Normalisation 

In this Appendix, we derive the scattering boundary conditions required above. 
The scattering theory is determined by a relation of form 

an(-^)l"+*lj|n+$|(fcr)+/3„(-z)-l"+*lj_|n+$|(fcr) ''^ (_z)V„(A:r) + ^^ (32) 

where the /„ are the scattering amplitudes and the integer order Bessel functions 
arise from the expansion of the incoming plane wave. All Bessel functions may 
be replaced by their asymptotic forms J^(x) ~ \J2/ {nx) cos(x — |)7r/2). The 
scattering normalisation is then determined by requiring the coefficients of e"*'^^ 
to match and gives + 1- This leads to the normalised spinor (|13]) given 
in the text. 



B Bounds on Scattering Lengths 

We derive here the range of allowed sector by sector. The scattering lengths 
R~ may be derived from these by -R^($) = -Rli_„(— Equation (|^ at zero 
kinetic energy (u;^ = im?) is 

'^^-f^^±^V + ^^U = 0. (33) 



r dr dr 



This factorises as 



(d^ n + l + air) \U_ n + air) \^^ (34) 
\dr r J \dr r J 

and so may be solved exactly to give two independent solutions 

^W(r) = r''expy\{r')/r'dr'^ (35) 

^(2)(r) = r"exp||J'^c/r'|^V'-'-'"exp|-2^' ^^rfr''|c/r^(36) 

The solution Fn{r) to (^ at zero kinetic energy is the solution with regular 
boundary conditions at r = 0. For n >0 this is clearly ip^^^r) whilst for n < 0, 
ip^"^^ is the appropriate solution. 
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Case (i): n > For r > a, ip^^\r) oc r"+*. Comparing with (|T9|) or using 
the fitting formula (|20|) , we see that = if n > — [$] (matching to r'""^*!) or 
= oo otherwise (matching to r"'""*"*'). 

Case (a): n < In general, ijj^'^\r) matches to a non-trivial linear combination 
of r^'""*"*' so the situation is more complex. The logarithmic derivative is given 

by 



F' 



n + $ f(a 



where /(r) is defined by 

fir) 

and /(r) by 



r ^" exp 



af{a) 

" a{r') 
r' 



dr' 



1 r 
r Jo 



f{r')dr' 



(37) 



(3^ 



(39) 



Clearly, /(r) and /(r) are positive and non-vanishing except at r = 0. We may 

\n + $1 



rewrite the fitting formula (20) as 



(i?+/a)2|"+*l = 1 
Thus if n > — [$], we obtain 

Rt = c 



2\n + $ 



+ 



F' 



(40) 



/(«) 



2\n + ^\f{a) 



l/(2|n+*|) 



(41) 



from which we can conclude the bound < R^ < a, given our observations about 
/(r) and f{r). On the other hand, if n < — [$] we find 



R' 



1 - 2|n + $ 



M 

/(«) 



l/(2|n+$|) 



(42) 



yielding the bound -oo < (i?+)2''^+*l < a2''^+*l. Note that R+ must be finite as 
a consequence of the non- vanishing of /(a). 

Furthermore, one can show that the above bounds are best possible in the 
sense that, for given n and any radius a, there exist magnetic flux distributions 
supported within radius a with any scattering length in the above allowed ranges. 
This may be proven by observing that a potential V{r) takes the form (|^) with 
a{r) smooth and obeying a{r) = Ar^ + O(r^) as r — > and a{r) = $ for r > a, 
if and only if the equation 



1 d d ^\ 

— r^-r + y ] u 
r dr dr / 







(43) 
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has a smooth solution u{r) which is non-vanishing in (0, oo) and obeys u{r) = 
r"(l + Ar^ + 0(r^)) as r ^ and u{r) oc r""*"* for r > a, whereupon we may 
identify 

n. (44) 



u 

air) = r — 
u 



Full details will appear in [Q. 



We remark that in |jT^ it is mistakenly concluded that in the limit as a — > 0, 
the upper spinor component is always either J\n+'i>\ or J-|„+$| and that (translated 
into our language) i?^ ^ if $ > in the critical sector. However, our arguments 
above and in |^ hold for arbitrarily small a and thus the range of scattering 
lengths allowed in the limit as a — is simply the appropriate limit of the 
range for finite a. Thus in the critical sector for $ > 0, any scattering length 
in the range — oo < (-R^)^'""'"*! < is allowed as the limit of scattering lengths 
of a sequence of flux tubes of decreasing radius. In we also prove rigorous 
statements about the convergence of the associated sequence of Hamiltonians. 
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T7 > n 


i?+ — 


-oo < < „2|n+*+l| 


-[<!>]< n < -1 


< it!+ < a 




n = -1 - [$] 


-oo < (i^+j^l'^+^l < a2|"+*l 


i?- = oo 


n < -2 - [$] 







Table 1(a): Allowed scattering lengths for $ > 0. 



n > -[$] 


i?+ = 


-oo < < a2|"+*+i| 


n = -1 - [$] 


i?+ = oo 


-oo < < a2|"+*+i| 


< n < -2 - [$] 


i?+ = oo 


< i?;; < a 


n < -1 


-oo < < a2|"+*l 





Table 1(b): Allowed scattering lengths for $ < 0. 
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3-2-10123 
Fig. 1(a) Maximum and minimum amplification factors. 




3-2-10123 
Fig. 1(b) Amplification assuming GUT scale scattering lengths. 

p($) ' 




3-2-10123 



Fig. 1(c) Amplification from the critical sector with GUT scale scattering 

lengths. 
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